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During the growth of metal thin films on dielectric substrates at a given deposition temperature X^, the film's 
morphology is conditioned by the magnitude and asymmetry of up- and downhill diffusion. Any severe change 
of this mechanism leads to a growth instability, which induces an alteration of the thin film morphology. In 
order to study this mechanism, ultra-thin Pt films were deposited via pulsed laser deposition (PLD) onto 
yttria-stabilized-zirconia single crystals at different deposition temperatures. The morphological evolution of 
Pt thin films has been investigated by means of scanning electron microscopy (SEM), atomic force microscopy 
(AFM) and standard image analysis techniques. The experimentally obtained morphologies are compared to 
simulated thin film structures resulting from a two-dimensional kinetic Monte Carlo (KMC) approach. Two 
main observations have been made: i) thin Pt films deposited onto Zr0 2 undergo a growth transition from 
two-dimensional to three-dimensional growth at > 573 K. The growth transition and related morphological 
changes are a function of the deposition temperature, ii) A critical cluster size of i* =4 in combination with an 
asymmetric Ehrlich-Schwoebel (ES) barrier favoring the uphill diffusion of atoms allows for a computational 
reproduction of the experimentally obtained film morphologies. 
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I. INTRODUCTION 

Pt thin films are commonly used as gate material on 
metal-insulator- metal (M-I-M) devices 1 , auto-catalysts 2 
and electrodes in micro-solid oxide fuel cells 3 5 . How- 
ever, thin metal films on dielectric substrates are ther- 
modynamically instable. The instability originates from 
the conflictive bonding characteristics of metals and in- 
sulators , which may be expressed by the interfacial en- 
ergy of metal/dielectric interface, and from the thin film's 
metastable configuration as a consequence of the chosen 
deposition method. This is especially true for deposition 
methods like PLD or sputtering, for which the kinetic en- 
ergy £?kin of the deposited atoms generally exceeds their 
thermal energy .Ether- For the deposited film, this results 
in a metastable configuration that tends to equilibrate, 
once subjected to temperature by annealing, joule's heat- 
ing or radiation. However, if it was possible to keep the 
deposition rate low, stable thin film structures artificially 
tailored for specific applications could be achieved in de- 
pendence of the deposition temperature T^. This requires 
the control of the instantaneous deposition flux (initial 
nuclei density), the average deposition flux (growth rate) 
and the deposition temperature (diffusion) . These pre- 
requisites are met for metal growth by pulsed laser depo- 



sition (PLD) 8 ' 9 , so called laser-MBE, especially in case of 
the large ratios between the pulse duration (« 10 — 25 ns) 
and the repetition rate (~ 1 — 100 Hz). This circumstance 
makes PLD ideal for surface-engineering 10 and for study- 
ing the fundamentals of metal-on-insulator growth. In 
order to gain a general understanding of the main prin- 
ciples of metal-on-insulator growth, Pt thin films on sin- 
gle crystalline yttria stabilized zirconia (Zr0 2 ) have been 
chosen as model metal/ceramic systems in the present 
study. The investigation will focus on how the deposi- 
tion temperature affects the operative growth mechanism 
which determines the final thin film morphology. In or- 
der to identify the dominating atomistic processes dur- 
ing growth the experiments have been performed along- 
side with kinetic Monte Carlo simulations. This paper is 
structured as follows: Section I introduced the basic con- 
cepts of metal-on-insulator growth and its importance for 
various applications. Sect. II deals with the experimen- 
tal framework. In Section III, detailed experimental and 
simulation results are presented and discussed. The fi- 
nal Sect. IV encompasses a summary of the findings and 
conclusions. 



II. EXPERIMENTAL 
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A. Sample Preparation 

In order to study the temperature dependence of 
the growth kinetics of metals on insulators, ultra-thin 
Pt films were deposited onto single-crystal yttria- 
stabilised-Zr0 2 (001) via PLD (F = 4.0(2) J/cm 2 , 
Pbase = 1 • 10" 6 mbar, p Ar = 2.7 • 10" 2 mbar). F 
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FIG. 1. Schematic illustration of the basic hopping processes 
used for the KMC model simulating the growth of Pt(lll) 
on Zr0 2 . (a) The growth rate is 1 atom per 2000 diffusion 
steps (exp. 1:4- 10 6 ) with a critical cluster size i* — 4 (b) 
Main processes A: group of atoms split, B: atom diffusion on 
Zr0 2 , C: atom going up a step, D: atom diffusion on Pt step 
(c) Temperature-dependent fall or jump height. 



denotes the laser fluence and p the pressure in the 
vacuum chamber. The substrates were precleaned using 
acetone and isopropanol and were thermally smoothed 
in a muffel furnace for 4 h at 1673 K before deposition. 
The Pt target (purity ~ 99.98 %) was ablated with a 
KrF + excimer laser with a wavelength of 248 nm, a 
pulse duration of 25 ns and a pulse repetition rate of 
10 Hz, leading to a ratio between pulse length and time 
of 1/4-10 6 . 

In order to achieve a measurement scheme of sufficient 
significance, Pt layers with different thicknesses (num- 
ber of shots n s = 25000 — 250000) were deposited at 
various deposition temperatures Td ranging from 473 to 
1073 K. It is noteworthy that the Pt/Zr0 2 is suitable 
as metal-on-insulator growth model system at elevated 
Td, because it is characterized by a chemically inert 
interface up to 1273 K 11 . 



B. Kinetic Monte Carlo Simulations 

A conventional two-dimensional kinetic Monte Carlo 
(KMC) approach has been chosen to reproduce the ex- 
perimentally observed 2D-3D growth transition. Thereby 
atoms are deposited on randomly picked sites on a pe- 
riodic grid with 3000 x 30 atomic positions. The model 
considers three main kinetic processes: (i) deposition of 
atoms with a sticking coefficient of a = 0.9, (ii) the site- 
specific surface diffusion of adatoms and (hi) the ascend- 
ing and descending of atoms from steps with different 
heights. 

The position-dependence of the diffusion rate is imple- 
mented using a variable energy barrier E^j, similar to 
the approach of Merrick et al. 12 . For a transition from a 
site i to a site j the diffusion rate is given by an Arrhenius 



TABLE I. Chosen barrier heights used for the various jump 
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law and reads as follows 
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(1) 



Thereby the attempt frequency z/q^ is chosen to be 
substrate-dependent and is for atoms on Pt z^o,Pt = 
1.2 -10 9 and ^o,Zr0 2 = 6.3 -10 9 for atoms on the Zr0 2 sub- 
strate respectively. The barrier height Ep th ^p t = 0.8 eV 
for Pt-Pt self-diffusion has been adapted from field-ion 
microscopy results , whereas the barrier height for the 
Pt diffusion on Zr0 2 has been set to £?zr0 2 ^Zr0 2 = 
1.0 eV. 

The dissociation of unstable clusters of atoms was taken 
into account by the introduction of a critical cluster size 
i* = 4. The energy barrier, which corresponds to the en- 
thalpy of the adatom formation, scales with the cluster 
size accounting for the increased nearest neighbor inter- 
action, see Fig 1 (a). Equivalent to the dissociation of 
unstable atom clusters, an Ehrlich-Schwoebel (ES) bar- 
rier is implemented in the KMC-model that accounts for 
an asymmetric probability of ascending and descending 
steps. A detailed overview of the different energy barriers 
chosen for the different jump processes is given in Tab. I. 
In addition to the ES barrier and in analogy to Leal 



et al., the step height h ste p t 



-14.2 + 0.03 • T d , 



that an atom can overcome by uphill diffusion scales lin- 
early with deposition temperature in the interval Td = 
[473, 1073] K, see Fig 1 (c). However the step height an 
atom can fall down is kept constant h step j. = 2 for all 
Td > 473 K. The same is true for internal jumps on Pt 
clusters. 



C. Sample Characterization 

The morphology of the samples was analyzed via high 
resolution AFM, using a Mobile S (Nanosurf) and a 
Topometrix 2000, and standard scanning electron mi- 
croscopy (SEM). The coverage, the distribution of nucle- 
ation centers and the island distribution were calculated 
from SEM or AFM images using standard particle anal- 
ysis algorithms. While the film thickness has been quan- 
tified by Rutherford backscattering spectrometry (RBS), 
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FIG. 2. Observed growth modes of Pt thin films deposited at different Td captured by AFM. (a) 2D layer-by-layer growth 
mode with a film thickness h — 16 nm (determined by RBS). (b) and (c) 3D kinetically frozen island growth characterized 
by a percolating network of elongated islands, (d) 3D island growth mode with isolated islands having a mean perimeter of 
120(50) nm. 



the orientation of the grown films has been determined 
using grazing incidence X-ray diffraction (GIXD). 
The simulated profiles and structures were compared to 
the films prepared by PLD using common statistical tools 
for surface description, i.e., the roughness and the height- 
height correlation function (HHCF). 
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For each output file of the simulation, the average rough- 
ness R a and root mean square roughness R rms were cal- 
culated according to Eq. (2) and Eq. (3), as well as the 



HHCF, 
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Thereby, d = Xi+i — X{ is the distance between neighbor- 
ing points in units of the grid. The determined HHCF is 
compared by a least square fitting to a Gaussian-HHCF 
(Eq. 5), which is characterized by the root mean devia- 
tion of heights a and the correlation length L c . 
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FIG. 3. Surface coverage G as function of the deposition 
temperature Td. The coverage decreases linearly for Td > 



573 K with a = 8.3(5) • 10" 4 K" 



III. RESULTS AND DISCUSSION 
A. Nucleation and Growth 

Pt thin films were grown on yttria-stabilized-Zr0 2 
single crystals at a pulse frequency of / p = 10 Hz 
over a range of different deposition temperatures, Td = 
[473, 1073] K and deposition times of 25000 - 250000 
shots. XRD measurements have shown that the Pt thin 
film preserves a (lll)/ cc / /(001) sc orientation relation to 
the Zr0 2 substrate for all Td. The overall growth rate 
is not affected by Td and is 1.1(5) ML (monolayers) per 
1000 laser pulses or 0.19(8) nm/min. For all deposited 
Pt thin films, the morphological transition as a function 
of the deposition temperature Td has been analyzed ex 
situ via AFM and SEM. 



1. Temperature domain 

In Fig. 2, the experimentally observed three main 
growth regimes of Pt on Zr0 2 is shown. For sufficiently 
low deposition temperatures Td < 473 K, the thin film 
grows layer-by-layer, see Fig. 2(a). The film thickness h, 
determined via RBS, is 16 nm and exceeds significantly 
the measured height profile of roughly 4 ML in Fig. 2(a). 
Therefore, the observed islands can be associated to a 
yet uncompleted layer of the thin film. 
Choosing Td > 573 K, the film changes its growth 
regime from a 2D to a 3D growth and undergoes a 
continuous-discontinuous transition. This is indicated 
by elongated clusters of grains and holes that are formed 
(Fig. 2(b)) which breaks the symmetry of the film, see 
Fig. 2(c). Increasing Td further results in an increase in 
the mean island size and, as a consequence, in a linearly 
decreasing coverage, see Fig. 3. In the last stage, the 
discontinuous film loses its percolation and an ordered 
array of faceted particles is formed as shown in Fig. 2(c). 



In the early stages of film growth, i.e. for Td < 573 K, 
the surface diffusion and surface self-diffusion of Pt is 
suppressed. This results in a decreased critical diffusion 
length l c (the mean distance between an adatom and a 
pre-existing nucleation site) and an enhanced density 
of nucleation sites due to the decreased mobility. It is 
noteworthy, that the coalescences' kinetics of different 
stable islands is decelerated as the capillary forces 
depend on the surface diffusion D s too . Therefore, the 
found layer- by- layer growth has to be understood as a 
result of kinetic processes rather than thermodynamic 
relaxation. Based on the assumption that the density 
of nucleation clusters p nuc is enhanced and the cluster 
coalescence is negligibly slow, the probability that atoms 
are deposited beyond the edge of steps is increased. In 
combination with preferential downhill funneling 15 ' 16 at 
these step edges, this leads to a layer-by-layer growth 
and a smooth film topography as shown in Fig. 2(a). 
Once the substrate temperature is shifted to higher 
temperatures, Td > 573 K, the surface diffusion is 
increased and the capillary forces gain more impact 
on the morphological evolution of the thin film. This 
onset of surface diffusion is featured in the observed 
2D-3D growth transition, which can be seen in Fig. 2(a) 
and Fig. 2(b). In the early stage of 3D growth, the 
formed film morphology can be understood as result 
of two competing mechanisms: the radial expansion 
of single clusters and the coalescences of two or more 
impinging clusters 17 . If the expansion kinetics of 
those clusters is much faster than the coalescence, the 
growth is kinetically frozen 18 and pores form within a 
percolating network of elongated grains, see Fig. 2(b). 
With increasing substrate temperature Td > 773 K, the 
porosity and the island size, with a mean perimeter 



of dj 



land 



120(50) nm, increase (Fig. 2(c)) until the 



percolation of the thin film is lost (Fig. 2(d)). These 
findings are in good accordance with the morphological 
evolution of epitaxial growth of Ag on Mica for different 
deposition temperatures 19 . The present morphologies 
and especially the decreasing coverage with increasing 
Td can well be explained under the assumption that 
adatoms can jump from the substrate on pre-existing 
clusters. In a first approximation, the thermal activation 
of the jump height is assumed to scale linearly with the 
film coverage (Fig. 3). It is noteworthy that the present 
morphological evolution at high Td is presumably due 
to the same capillary forces, that cause the continuous- 
discontinuous transition of Pt on Zr0 2 during thin film 
agglomeration 20 . 



2. Time domain 

In order to compare the obtained growth morphologies 
with findings on metal-on-insulator growth via PLD in 
literature, a set of films were deposited at Td = 773 K for 
three different deposition times. The resulting morphol- 
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FIG. 4. SEM images of Pt film growth at T d — 773 K for different times (number of shots), (a) Nucleation of islands on a 
vicinal terraced substrate with enhanced nucleation at the terrace steps. The mean terrace width is dkink = 245(8) nm. (b) 
and (c) 3D kinetically frozen island growth, characterized by a percolating network of elongated islands. 



ogy evolution as function of deposition time, i.e. number 
of shots, is depicted in Fig. 4. In Fig. 4(a), the initial 
stage of film growth is shown. Isolated Pt nuclei with 
a mean next-neighbor distance of l e ff = 32(1) nm have 
been formed on the terraces and along the step edges 
of the vicinal Zr0 2 substrate. The density of nuclei per 
unit area was determined to be p nuc = 609 fim~ 2 . The 
enhanced growth at the step edges is in good accordance 
with results from literature on the epitaxial growth of 
Pt (111) on (0001) sapphire 21 ' 22 . In these cases, it has 
been found that the sapphire surface steps act as nucle- 
ation sites and cause rotational twinning in the Pt-film. 
With increasing deposition time, the nuclei grow and co- 
alesce with each other, see Fig. 4(b). Because of the 
different velocities of the cluster growth and the coales- 
cences of clusters the resulting film is disrupted by holes, 
i.e. the growth is kinetically frozen . Further deposi- 
tion as shown in Fig. 4(c), increases both the feature size 
of impinging grains and the hole size. In addition, new 
nucleation clusters form inside the cavities of the dis- 
rupted film. The found growth morphologies and their 
progression with increasing deposition time are in good 
agreement with the experimental findings on metal-on- 
insulator growth of Ag on Mica 18 ' 23 ' 24 . 



B. Experiment vs. Simulation 

In order to substantiate the previous findings, the sig- 
nificance of an increasing uphill diffusion with increas- 
ing deposition temperature Td regarding the morpho- 
logical progression of Pt on Zr0 2 has to be confirmed. 
In Fig. 5, the morphological regimes during thin film 
growth, simulated by the two-dimensional KMC model, 
are shown. The chosen approach reproduces the experi- 
mentally observed transition from layer-by -layer growth 
to 3D mound formation at Td ~ 573 K. The continuous- 
discontinuous transition of the thin film for Td > 673 K 
has been reproduced as well. Thus, it can be concluded 
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FIG. 5. Surface morphologies obtained from KMC simula- 
tions after 20000 growth steps for 6 different deposition tem- 
peratures Td. The growth conditions as well as the obtained 
surface morphologies scale with altered substrate temperature 
T d . 



that the experimentally observed growth transitions can 
be treated in terms of a thermally activated step height 
dependent ES barrier, whereby the jump height of a 
diffusing adatom increases linearly with Td. In order 
to quantitatively validate the morphological resemblance 
between the experimental and simulation results, the 
HHCF was chosen measuring the scaling properties and 
roughness for all obtained thin films. The HHCF cap- 
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FIG. 6. (a) Plot of the root mean square roughness R rms VS. 
Td for the experimentally grown and the simulated Pt thin 
film morphologies, (b) Plot of the correlation length L c as 
function of T^, comparing the experimental with KMC simu- 
lation results. The gray shaded area denotes the 2D growth 
regime 



tures two main surface characteristics, the high frequency 
fluctuations perpendicular to the surface represented by 
the root mean deviation of heights a and the correla- 
tion length L c . The latter denotes the critical distance, 
within characteristics of the surface can be expected to 
be correlated !5 ' 26 . The measured RMS-roughness for 
the experimental and the simulated surfaces is shown 
in Fig. 6(a). The step-like profile coincides with the 
observed continuous-discontinuous transition driven by 
the enhanced capillary forces. The RMS-roughness for 
the experimental and the simulated surfaces follows the 
same trend. The observed difference by a factor of four 
between the experimental data (deposition on a 2D sur- 
face) and simulated data (deposition on a ID line) is 
attributed to dimensional effects. A similar conclusion 
can be drawn regarding the evolution of L c as a function 
of the deposition temperature T^: L c is in the order of 
the film's grainsize for low and of the island size for 



high Td. The surface characteristics of both experimental 
and simulation results are in good accordance, so that it 
can be concluded that the main growth mechanisms have 
been found for the present metal- insulator couple. 



IV. CONCLUSION 

In essence, it has been shown that Pt thin films de- 
posited by PLD undergo a 2D-3D growth transition as 
function of the deposition temperature X^. The growth 
kinetics and morphological evolution depend strongly on 
the onset of the surface self-diffusion and its related cap- 
illary forces. Three main growth regimes have been iden- 
tified: i) layer-by-layer growth governed by downward 
funneling due to the absence of capillary forces, ii) ki- 
netically frozen island growth, which is attributed to 
the onset of surface diffusion that competes with nuclei 
growth and coalescence and iii) pure 3D island growth, 
which is dominated by uphill diffusion due to capillary 
forces. The choice of the appropriate deposition tempera- 
ture thus allows for an application-specific tailoring of the 
thin film's structure and hence its properties. In particu- 
lar, these morphological transitions are in excellent agree- 
ment with the results of temperature-dependent growth 
experiments of Ag on Mica by Baski and Fuchs . Addi- 
tionally, the temperature-dependent growth kinetics and 
morphological evolution of ultra-thin Pt films have been 
investigated numerically using Kinetic Monte Carlo sim- 
ulations. The simulations were performed using a critical 
cluster size of i* =4 and an asymmetric, temperature- 
dependent ES barrier facilitating the upward diffusion 
of Pt atoms on pre-existing terraces. This allowed for 
a qualitative reproduction of the experimentally deter- 
mined surface roughness and height height correlation 
length, which confirms the pre-dominant impact of cap- 
illary forces on the thin film growth at elevated depo- 
sition temperatures. It has been shown that a detailed 
knowledge of the interacting physical growth mechanism 
is essential to control the engineering of thin films at el- 
evated deposition temperature. Within this manuscript 
the capillary forces, which also control solid state dewet- 
ting phenomena 20 , have been identified as the most cru- 
cial. 
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